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We consider an Ising model on a square lattice with ferromagnetic spin-spin interactions spanning 
beyond nearest neighbors. Starting from initial states with a single unbounded interface separating 
ordered phases, we investigate the evolution of the interface subject to zero-temperature spin-flip 
dynamics. We consider an interface which is initially (i) the boundary of the quadrant, or (ii) the 
boundary of a semi-infinite stripe. In the former case the interface recedes from its original location 
in a self-similar diffusive manner. After a re-scaling by \/t, the shape of the interface becomes more 
and more deterministic; we determine this limiting shape analytically and verify our predictions 
numerically. The semi-infinite stripe acquires a stationary shape resembling a finger, and this finger 
translates along its axis. We compute the limiting shape and the velocity of the Ising finger. 
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I. INTRODUCTION 

The concept of universality plays a crucial role in un- 
derstanding of myriads of many-particle systems. Uni- 
versality is the claim that details are irrelevant, so that 
an interesting behavior of a given system can be extracted 
from its simplified versions. In the context of Ising ferro- 
magnets, for instance, universality implies that the crit- 
ical behavior depends on the spatial dimensionality, but 
not on the details of the underlying lattice. Universality 
is also embodied in another even more 'obvious' claim, 
viz. that the toy Ising model with nearest-neighbor (NN) 
interactions faithfully represents the behavior of the gen- 
eral model with ferromagnetic interactions spanning be- 
yond nearest neighbors. Here we probe the validity of 
this latter assertion. 

Equilibrium behaviors at the critical temperature 
agree with universality when the interaction strength 
quickly decreases with separation between the spins pQ. 
Dynamical behaviors are critical [2 -4 J , as it is manifested 
e.g. by algebraic time dependencies both at the critical 
temperature and at sub-critical temperatures. Here we 
consider the latter, more specifically a zero-temperature 
spin-flip dynamics. We focus on the two-dimensional fer- 
romagnetic Ising model. On the square lattice, the small- 
est perturbation of the Ising model with NN interactions 
accounts for next-nearest-neighbor (NNN) interactions. 
The Hamiltonian is 



J£ = - J x s i s i - J2 
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where s — ±1 are Ising spins. Couplings in Eq. ([!]) 
are presumed to be ferromagnetic, J\ > and J2 > 0, 
but otherwise arbitrary. The first sum in Eq. is 
taken over NN pairs, that is, over sites i = (11,12) 
and j = (jijjz) such that i — j is equal to (±1,0) or 
(0, ±1). The second sum in the Hamiltonian is taken 
over NNN pairs of spins, more precisely over diagonal 
neighbors, that is, i — j attains one of the four values: 
(1,1), (1,-1), (-1,1), (-1,-1). _ 

At zero temperature, energy raising spin flips are never 



performed. A zero-temperature spin-flip dynamics im- 
plied by the Hamiltonian quantitatively differs from 
the standard dynamics. We are interested in qualitative 
changes, however. Surprisingly, such changes do exist, al- 
though they are rather subtle. For instance, starting from 
above the critical temperature and suddenly quenching 
to zero temperature, the standard Ising model on the 
square lattice can get trapped in a state with vertical or 
horizontal stripes [SH7]. For the Ising model character- 
ized by the Hamiltonian 0, stripes going in the (1, 1) 
and (1,-1) directions (the boundaries of such stripes are 
perfect ladders with steps of the size of the lattice spac- 
ings) are also possible [5]. Increasing the span of inter- 
actions results in positive (in the thermodynamic limit) 
probabilities of observing more exotic stripe states [5]. 

In this work we focus on simple deterministic initial 
settings, specifically on the behavior of a single interface 
separating ordered phases. One example is the evolution 
of a quadrant of minority phase (Fig. [lj. Both for the 
standard Ising model and for the model with Hamiltonian 
the major features of the dynamics are the same, e.g. 
the typical size of the 'melted' region grows diffusively, 
i.e., as \ft. After re-scaling by y/i, the boundary of the 
melted region becomes more and more deterministic and 
it approaches a limiting shape in the large time limit. 

We will show that the limiting shape arising for the 
Ising model with Hamiltonian differs from the lim- 
iting shape corresponding to the standard Ising model. 
Increasing the interaction span leads to a zoo of limiting 
shapes. To define the interaction range we endow the 
lattice (the square lattice in our case) with a metric and 
postulate that spins separated by distances exceeding the 
interaction range don't interact. The Manhattan metric 
|i — j| = l*i — ill + N2 — J2I is perhaps the simplest. We 
shall use it and we consider a class of Hamiltonians 
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Coupling constants are always assumed to be ferromag- 
netic, J n > for all n = 1, . . . , k, and rapidly decreasing 
with n. Once such requirements [([3| and (49) below] are 
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FIG. 1: The interface is the boundary of a quadrant in the 
initial state (left); the interface encloses an area St at a later 
time t. 



obeyed, the precise values of the coupling constants are 
irrelevant and only the interaction range k plays a role. 

The rest of this paper is organized as follows. We study 
in detail the Ising model which, in addition to the stan- 
dard NN interactions, has NNN interactions. It is easier 
to think about the Hamiltonian ([!]) than the Hamilto- 
nian ^ with k — 2, although there is no difference, e.g. 
in the case of a single interface shown in Fig. 
Hamiltonians lead to the same dynamics. In Sec. 
specific example of the 'melting' of the quadrant of mi- 
nority phase, we map the interface onto a lattice gas. We 
utilize this mapping in Sec. |III| We first examine the evo- 
lution of the density of the lattice gas, and then translate 
this knowledge into a calculation of the limiting shape of 
the interface. In Sec. |IV| we study how an external mag- 
netic field, which disfavors the minority phase occupying 
the quadrant, changes the limiting shape. The limiting 
shape of the finger is described in Sec.|V] All the limiting 
shapes in sections |III|[V| are derived in the realm of the 

we 



Ising ferromagnet with NNN interactions. In Sec. VI 



analyze larger interaction ranges, viz. we consider a class 
of Hamiltonians ( p|) w ith k > 3. Concluding remarks are 
presented in Sec. |VII| 



II. MAPPING OF AN INTERFACE ONTO A 
LATTICE GAS 

The zero-temperature dynamics allows spin flips only 
on interfaces separating ordered phases. This dynam- 
ics tends to preserve the integrity of a single interface. 
Thus, the two-dimensional problems essentially reduce 
to one-dimensional problems with all the 'action' hap- 
pening along the interface. The above assertions are not 
manifestly true. Consider a spin surrounded by four sim- 
ilarly aligned spins in an infinite sea of oppositely aligned 
spins. The interactions of the central spin with its four 
nearest neighbors must outweigh its interactions with the 
rest of the misaligned spins to assure that the central spin 
would not flip. This happens when 




n>2 



(3) 



FIG. 2: An illustrative interface rotated by 7r/4 and the cor- 
responding particle configuration. A spin-flip event is shown 
together with the corresponding hop of the particle in the 
lattice gas. 



An assumption that the coupling constants rapidly de- 
crease with n presumes the validity of the constraint 
([3]) which obviously guarantees that any spin inside a 
droplet of aligned spins cannot flip. The integrity of a 
single interface is not satisfied for finite interfaces, and 
even for semi-infinite interfaces bordering regions of fi- 
nite width (see Sec. [V] for more details); this already oc- 
curs for the Ising model with zero-temperature dynamics 
associated with the standard NN interactions. If, how- 
ever, the size of the interface greatly exceeds the lattice 
spacing, small (and quickly disappearing) droplets very 
rarely calve from the interface. In our first example of 
the melting of a quadrant, the interface manifestly keeps 
its integrity throughout the evolution. 

The dynamics greatly simplifies when we can map the 
interface onto a one-dimensional lattice gas. Such a map- 
ping is applicable throughout the evolution of a quadrant 
of minority phase (Fig. [IJ. Zero-temperature spin- flip 
dynamics leads to the melting of the quadrant, and the 
boundary of the melting quadrant remains an unbounded 
non-self-intersecting connected interface. The dynamics 
is stochastic and the interface is also stochastic, yet in 
the long time limit relative fluctuations vanish and the 
notion of the limiting shape becomes well-defined. The 
goal is to determine the limiting shape of the interface. 

The mapping of the interface dynamics onto a one- 
dimensional lattice gas is obtained by rotating the in- 
terface counter-clockwise by angle 7r/4 around the ori- 
gin and projecting the boundary onto the horizontal line 
(see Fig. [2]). We now identify bonds of the interface with 
sites on a one-dimensional lattice and we put a particle 
on a site (leave a site empty) if the corresponding bond 
on the interface goes along the co-diagonal (diagonal). 
There could be at most one particle per site (exclusion 
property). The hopping rates are affected by nearest- 
and next-nearest neighbors on the one-dimensional lat- 
tice in the case where the spin interactions are described 
by the Hamiltonian (|2) with k = 2; for an arbitrary fc, 
particles separated by distance < k from the particle that 
attempts a hop influence the dynamics. 

More precisely, in the particle mapping for the spin 
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Hamiltonian ^ with k = 2, particles can hop to neigh- 
boring empty sites (symmetric hopping), and such hops 
are carried out at random according to the following 
rules: 

1. If a hop would not increase the number of NN pairs 
of particles, it is always performed. 

2. If a hop would increase the number of NN pairs of 
particles, it is never performed. 

The above hopping rules directly follow from the zero- 
temperature spin-flip dynamics compatible with Hamil- 
tonian ([TJ or Hamiltonian ([2| with k = 2. The major 
property of the zero-temperature dynamics is that a spin 
flip is forbidden if it would lead to an increase in energy. 
There is a bit of arbitrariness in hopping rates, e.g., hops 
which decrease the number of NN pairs of particles may 
occur at a different rate than hops which do not change 
the number of NN pairs of particles. In the long time 
limit, however, the fraction of hops of the first kind be- 
comes negligible, and therefore it suffices to study the 
simplest version in which these rates are equal. We shall 
always set this rate to unity. 

We now illustrate the hopping rules by the early-time 
dynamics. We start with the unperturbed quadrant, or 
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in the lattice gas representation. The spin at the corner 
of the quadrant is the only flippable spin. Its flipping is 
equivalent to the hopping 



> o o o o o , 



• •••o»oooo. 



which is the only move compatible with the exclusion 
property. For the standard Ising model with NN spin- 
spin interactions, there are three possibilities for the next 
move: 



••••o»oooo 



• •••ooooo 

ItlOltOOOO 
• •IIOOIOOO 
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which all occur with equal rates. This illustrates the gen- 
eral assertion that the underlying lattice gas is a symmet- 
ric simple exclusion process [U IMTT] . For the model ([I]) , 
however, there are only two possible moves: 



• •••0*0000 



• **o**oooo 

• •••00*0 00 



(6) 



Returning back to the original quadrant is impossible 
since such a flip would raise the energy; according to the 
lattice gas description, such a move would have created 
an additional NN pair of particles, and therefore it is 
never performed. 

Therefore our lattice gas (we shall call it a repulsion 
process) is governed by the zero-temperature dynamics 
associated with the Hamiltonian 



niUi+i 



(7) 



Here rii = 1 if site i is occupied, and n, = if it is empty. 
The dynamics is based on the exchange, O rij+i. This 
exchange is forbidden if it would lead to an increase in 
the number of NN pairs of particles, and thereby to an 
increase in energy. 



III. MELTING OF THE QUADRANT 

Here we consider the Ising ferromagnet with NN and 
NNN spin interactions. We apply the above lattice gas 
representation of the boundary of the corner to the de- 
termination of its limiting shape. The first step is to 
calculate the average density of the lattice gas, which is 
then used (through the inverse mapping from the lattice 
gas to the interface) to extract the limiting shape. 



A. Average Density 

The density rij(t) is a strongly fluctuating quantity. 
On the hydrodynamic level we are interested in distances 
greatly exceeding the lattice spacing. In this situation we 
can talk about the average density p(z,t). This quantity 
obeys a macroscopic diffusion equation 



dp 

at 



d_ 

dz 



D(P) 



dp 



The initial condition Q becomes 



p(z,t = 0) 



1 z < 
z > 



(8) 



(9) 



on the macroscopic (hydrodynamic) level. 

The diffusion description Q has been justified for nu- 
merous lattice gas models O |13j , but the diffusion co- 
efficient D(p) has been computed in very few models, 
mostly for lattice gases with density-independent diffu- 
sion coefficient. For the repulsion process, the diffusion 
coefficient depends on the density. This already follows 
from a few exact values: 



D(P) 




(10) 



The predictions of Eq. ( 10 1 can be understood as follows. 



In the small density limit, we can consider one particle 
in a vacuum. This particle diffuses with D = 1 (recall 
that we have set the hopping rates to unity). Similarly 
when the density is very close to the maximal density 
p = 1, we can consider a single vacancy in an otherwise 
fully packed lattice. This vacancy diffuses with D = 1. 
Thus D(p = 0) = D(p = 1) = 1. Consider now a small 
perturbation of the p — \ state. There is just one half- 
filled equilibrium configuration, ...•o»o»o»o.... A 
perturbation of the p = \ state is the configuration 



•o»o»o»»o»o»o»o 
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with one doublet. In this configuration both particles 
forming the doublet can hop. A single hopping event, 
say to the right, leads to 



we recast ( 16 ) into 



o • • o • o 



o • o • • o 



(11) 



Thus the doublet effectively hops with the same rate 1 
as the particles, and since the span of the doublet hops 
is twice larger (two lattice sites), the diffusion coefficient 
is 4 times larger. These arguments explain (10 1. 



The diffusion coefficient of an interacting dense lattice 
gas is usually impossible to determine analytically. For- 
tunately, for the repulsion process with Hamiltonian (17j , 
the diffusion coefficient has been recently computed [14] 



D(p) 



P~ 2 



~ 2 0<p<i 
I <P< 1 



(12) 



The derivation [14] is lengthy. (One must understand 
the steady states of the repulsion process and then use 
this understanding to compute the free energy, the com- 
pressibility, and employ the Green-Kubo formula |12j to 
determine the diffusion coefficient.) Equation ( 12 ) agrees 



with ( 10 ) derived above, as well as with the mirror sym- 
metry D(p) = D(l — p). 

We now return to the boundary- value problem ([8])-([9]). 
The self-similar nature of the problem (viz. the lack of 
the characteristic spatial scale) tells us that we can seek 
a solution in a scaling form 



p(z,t) = /(0, c 



z 
It 



(13) 



The mirror symmetry allows us to limit ourselves to the 
p < J region corresponding to the half-line z > (ini- 
tially the vacuum state). Inserting (13 1 into ^ and us- 



ing Eq. (12) for the diffusion coefficient, we arrive at the 



ordinary differential equation for the scaled profile 



d 



(1-/) 



-2# 



2C|=0 
d( 



while the boundary conditions are 

f(0) = L /(oo) = 



(14) 



(15) 



B. Limiting Shape 

The variables x and y characterizing the interface in 
the original reference frame can be determined from 



Writing 



y = I dzp(z,t) 

lx-y 



x y 
U \/4i 



(16) 



(17) 



V 



dC/(0 



(18) 



t-ri 



which (implicitly) determines the scaled limiting shape 
of the interface. We solved the boundary-value prob- 
lem ( 14)-([l5"| numerically. The scaled limiting shape is 
plotted on Fig. [3] We also performed numerical simula- 
tions of the repulsion process, and we used the discrete 
version of Eq. ( |18| to obtain the corresponding limiting 
shape. We ran 100 realizations of the dynamics, each 
until time t — 2 24 , and we averaged the interface profile 
over all realizations. We do not plot simulation results 
since the discrepancy between the simulated growth pro- 
cess and the numerical solution of the boundary-value 
problem (14)-(15) is undetectable by eye. 




FIG. 3: The limiting shape of the melting quadrant corre- 
sponding to the Ising model with NNN spin-spin interactions 
(k — 2). The limiting shape arising in the standard Ising 
model (k — 1) is shown for comparison. 

We now discuss two basic features characterizing the 
interface — the intersection of the interface with the di- 
agonal and the area under the interface. 



1. Distance from the origin 

The diagonal, that is the ray in the (1,1) direction, 
crosses the interface at 

/>oo 

Cdiag = ?7diag = / dC/(C) (19) 
JO 



Integrating (14) from £ = to £ = oo, one finds 
dC/(0 = -2/'(0) 



(20) 
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Combining (17), (19 1, and (20 1, we arrive at 



Inner and Outer Corners 



with 



^diag 2/diag C V^j 



C = -4/'(0) = 0.8655. 



(21) 



(22) 



More precisely, Eq. (21) is the average position of the 



intersection of the interface with the diagonal. The same 
behavior (21) arises for the standard Ising ferromagnet 



with NN interactions, but in that case the amplitude is 
C = tt- 1 ' 2 = 0.564189583.... 



2. Average area 
The average area is equal to 
(St) 

It is convenient to re-write / as a double integral 



ydx = tt = itl (23) 

ii Jo 



d^drj = 2 



d£ drj 



£>rj>0 



Changing variables, (£,?y) — > (u = £ — ??,??), and noting 
that the Jacobian is equal to unity, gfe^j = 1, we get 



d£dr) = 2 r/du (24) 

£>J7>0 Jo 



Plugging (18) into (24) we get 



2 du 
la 



dC/(0 



2 / dCC/(C) 
lo 



(25) 



We now multiply Eq. ( 14 ) by £ and integrate from C, = 



to Q = oo. Using integration by parts we obtain 



2 / dCC/(C) 



1 

2 jo 
1 rV2 



daw)- 2 ! 



2 7 (i - /) 2 : 

Thus J = | and therefore the average area is 



(26) 



Interestingly, the average area in the case of the standard 
Ising ferromagnet, (St) = t, is twice smaller. 

In Fig. [3] we show the limiting shapes for the stan- 
dard Ising ferromagnet and for the Ising ferromagnet 
with NNN spin interactions. We re-scale the plots in 
such way that in the (£, rj) variables both areas are equal 
to unity: / °° rj d£ = 1. Hence instead of (17) the scaling 
variables are (£, 77) = i _1 / 2 (a;, y) in the case of the stan- 
dard Ising ferromagnet and (£, 77) = (2t)~ 1 ^ 2 (x, y) for the 
Ising ferromagnet with NNN spin interactions. 



Initially there is one inner corner and no outer corners 
(left part of Fig. [I]). Generally we denote by N + (t) [resp. 
N_(t)] the average total number of inner [resp. outer] 
corners. The conservation law 



N + (t)-N_(t)=l 



(27) 



has an elementary topological origin and it is always valid 
in two dimensions. The average numbers of inner and 
outer corners grow as y/i. This in conjunction with the 
conservation law (27) leads to conjectural behaviors 



N + = A 2 yft + B 2 
iV_ = A2V1 + B2 



1 



The leading asymptotic in ( 28a H 28b) is 



A 2 = 2C 



(28a) 
(28b) 



(29) 



with C given by (22). Indeed, let's start at the diagonal 
point and go to x — > 00. The interface will make down- 
ward steps of length one; more precisely, almost all steps 
will have length one in the long time limit. The reason is 
easy to see in the realm of the lattice gas representation. 
The gas is in local equilibrium, and equilibrium at a den- 
sity p < \ has a simple structure [14] : Adjacent particles 
are always separated by at least one vacancy. There- 
fore the total number of inner corners below the diagonal 
point is equal to ydiag in the leading order. Similarly, the 
total number of inner corners above the diagonal point is 
equal to Xdiag- Hence N + = x^mg + j/diag in the leading 
order, which implies (29). 



We do not have a rigorous justification of (28a)-(28b), 



although our simulations support these expansions. It 
would be interesting to determine B2 analytically. 



IV. INFLUENCE OF A MAGNETIC FIELD 

Here we consider the quadrant in an external magnetic 
field. We assume that the magnetic field favors the ma- 
jority phase (otherwise the quadrant does not evolve). 
For the Ising model with NN and NNN spin interactions, 
the Hamiltonian reads 



5£ 



-Jl ^ S i S J- J 2 s i s 3~ h Yl 



Si 



(30) 



|i-j|=l 



|i-j|=2 



with positive magnetic field, h > 0, if the majority phase 
is the plus phase (as in Fig. [I]). At zero temperature, the 
strength of the magnetic field is irrelevant as long as it 
is smaller than a certain threshold so that only spins at 
the inner corners can (in principle) flip. In this situation, 
in every move one minority spin flips. In the lattice gas 
representation, the dynamics is the biased version of the 
dynamics studied previously, viz. the hopping rules are 



() 



identical up to the additional constraint — only hopping 
to the right is allowed. 

The hydrodynamic description of the biased repulsion 
process is based on the continuity equation 



dp dJ_ 
di + ~dz~ 







The average current J(p) has been computed [14] 



m = 



' P(l-2p) 
l-n 
(l-p)(2p- 



< p < \ 
5 <P< 1 



(31) 



(32) 



The current is maximal when p* — A= and p* = 1 — 

with J(p*) = J(p*) = Jmax = 3 - 2V2. The relation 
between the maxima, p* + p* = 1, reflects the more gen- 
eral mirror symmetry, J(p) = J(l — p). Note that for the 
standard Ising model in a magnetic field the correspond- 
ing lattice gas is an asymmetric simple exclusion process, 
as was first discovered by Rost [IB] . 

The derivation of ( 32 ) is simple since the steady states 



of the biased repulsion process are known (they are the 
same as in the symmetric case) |14j . The expression 
(32 1 for the current, sometimes only in the high-density 
appears in a number of models 



regime p > | , 

We solve Eq. (31 1 by noting that the density admits a 
scaling form 



p(z,t) = R(Z), 



z= z - 

t 



(33) 



Plugging this scaling form into (31) and using Eq. (32) 



for the current we obtain a rarefaction wave 



R(Z) 




Z < -1 
-1 < z < 
< z < 1 
z > 1 



(34) 



when the initial condition is given by ^ 
lie 




The density profile of the rarefaction wave ( 34 ) is dis- 
with densities p(— 0) = p* - 



continuous at x 



and p(+0) 



s/2 



p* = 1 — y| , precisely the values where the 
current attains its maximum. The average total num- 
ber of particles that penetrates into the region x > is 
therefore N(t) — J max t, where J max ^3 — 2\/2 



The interface is determined from (16), which becomes 



Y = I dZR(Z) 

'x-Y 



in the scaled variables 



X = 



t 



Y = 



t 



(35) 



(36) 



Combining (|35j) with ( 34 1, we find that the limiting shape 



of the interface (see Fig. P4J) is given by 



Y(X) = 



1 - 2VX 0<X< J max 

(l-2X + X 2 )/4: J max <X<l 



(37) 




FIG. 4: Limiting shapes of the melting quadrant in a mag- 
netic field. Upper curve: The limiting shape y/X + y/Y = 1 
corresponding to the standard Ising model, see |15| . Middle 
curve: The limiting shape ( |37[ ) corresponding to k — 2. Lower 
curve: The limiting shape Eq. (62 I corresponding to k = 3. 



Alternatively, one can directly derive, see Appendix \K\ 
an equation describing the evolution of the interface in 
the long time limit (when we can ignore fluctuations): 



Vt = 



1 + fe)- 1 



> y x > 

-1 > Vx 



(38) 



where y t = -gf and y x — g| . The structure of ( 38 ) sug- 
gests to seek a solution in the scaling form Y = Y(X). 
Plugging this form into (38), we obtain a differential 



equation for Y{X) whose solution is given by Eq. (37) 



V. ISING FINGER 

We now turn to the finger geometry (Fig. [5]). We as- 
sume that the minority phase initially occupies the semi- 
infinite region: y > and \x\ < L. In the interesting 
long time regime, more precisely when t 3> L 2 , the two 
corners of the initial finger interact and the finger relaxes 
to a limiting shape moving along the axis of the finger 
(in the y > direction) with some constant velocity. In a 
reference frame moving with the finger, the interface y(x) 
is stationary. We shall determine this stationary limit- 
ing shape by employing a continuum description which 
is valid when the width of the finger greatly exceeds the 
lattice spacing: L>1. 

During the evolution of the finger one rarely observes 
a fission of tiny domains of minority phase from the fin- 
ger. Apart from pathological initial conditions, e.g. those 
which contain 'tendrils' of width one (i.e. equal to the 
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FIG. 5: Schematic illustration of the evolution of a semi- 
infinite strip (a rectangular finger of the minority minus phase 
surrounded by the majority plus phase). On the right-hand 
side, the flip of the lowest minority spin, the tip spin in this 
example, is an irreversible process that causes the minimum 
height of the finger to advance by one. 





FIG. 6: An illustration of a rare emergence of a short-living 
tiny closed interface during the evolution of the finger. 



lattice spacing) or strips of width one, or rare separations 
of tiny drops (Fig. [6|, such break ups play a significant 
role only for narrow strips, L = 0(1). In the interesting 
situation of L 3> 1 we can ignore fission events and fluc- 
tuations and focus on the deterministic limiting shape. 

Due to symmetry we can limit ourselves to the region 
< x < L and y > 0. The governing equation for the 
interface shape y{x,t) is [see Appendix [A] 



Ut — Vxx J 



Vxx 

vl 



Vi 



< y x < l 

1 < Vx 



(39a) 
(39b) 



Consider an upward moving reference system in which 
the interface is stationary. In the bottom part of the 
interface, < x < V with yet unknown V , we should 
Therefore 



use (39a 



Vxx, y(0)=y x (0)=0, y x (L') = l (40a) 



In the top part, L' < x < L, we employ (39b I, so we have 
y x {L') = l, y{L)=y x {L) = ^ (40b) 



Vxx 

vl 



We must solve (40a) and (40b), match the solutions, and 



self-consistently determine the intermediate point L' and 
the velocity v of the finger. 



Transforming the variables 



X = 



L 



r = f, v = vL, e= L T 



and writing F = 4£ we recast (40a)-(40bl into 



V = 



dF 
dX 



F(0) = 0, F{1) = 1 



V = F~ 



dF 
dX 



F(£)=l, F(l) 



Integrating (42a)-(42bl we obtain 



F = VX, Vi = 1 
F = 1 



1 - V{X ~ i) 



(41) 

(42a) 
(42b) 

(43a) 
(43b) 



Combining (43b) with F(l) = oo and VI — 1 we get 

V = 2, 1=1 (44) 



Thus the velocity of the finger is 

2 

v = — 
L 



(45) 



For the Ising ferromagnet with NN spin-spin interactions, 
the velocity v = 1/L (see jTTj ) is twice slower. This two- 
fold acceleration is analogous to the twice faster growth 
of the area, Eq. (|26|, in the corner problem. 



Using Eqs. (|43a|)-(|43b|) and (|44|), we integrate F = f£ 



and find the shape of the finger(see Fig. [7]) 



Y 



X 2 



i-|ln(2-2|X| 



i < II|< l 



(46) 



VI. LONG-RANGED INTERACTIONS 

Here we study Ising ferromagnets with longer interac- 
tion range. Specifically, we consider the class of Hamil- 
tonians ([2| with k > 3. When k = 3, the corresponding 
lattice gas has NN and NNN repulsive interactions 



oo 

#3 = h X! n i n i+l 



oo 

- ^3 ^2 TliTli+Z 



(47) 



In the general case, 

oo oo 

H k = -h ^2 riim+i + ... + J k ^2 n t n t+k-i (48) 

i— — oo i— — oo 

Note that the interaction span for the repulsion process 
(481 is equal to k — 1, while in the underlying Ising fer- 



romagnet the interaction span between the spins is k. 

The correspondence between the Hamiltonians ^ of 
the two-dimensional Ising ferromagnets and the Hamilto- 
nians Q, (|47| and (|48]) of one-dimensional lattice gases 



8 



2.5 



2.0 



1.5 



Y 



1.0 



0.5 



1 I I I 


1 

— k = 


1 ! 




k = 


2 ; 


_\ \ 

\ \ 

\ \ 
_ \ \ 
\ \ 

\\ 
\\ 




/ / 
/ / 

i 

1/ 

u 


vv. 
\\ 
\\ 




A 
// 
// 


1 


- — i — 





-0.5 





X 



0.5 



FIG. 7: The limiting shape, Eq. (|46|, of the finger. The limit- 
ing shape Y — —\X\ — ln(l— \X\) (see [11]) of the finger in the 
case of the standard Ising dynamics is shown for comparison 
(dotted blue curve). 



is not straightforward, e.g. the Hamiltonians ([2| involve 
the interactions in the bulk (away from the interface) 
which are omitted. What is important is that the non- 
conservative spin-flip dynamics caused by ^ is the same, 
in terms of the interface, as the conservative hopping dy- 
namics implied by the Hamiltonian Q for k = 2, (47) 
for k = 3, and (48 1 for general k. 



Narratively, the dynamics corresponding to ( 47 ) is the 
following: 

1 . If a hop would decrease the number of NN pairs of 
particles, it is always performed. 



the right is allowed), the current J(p) is given by 

< p < | 



p(l-3p) 

l-2p 
(l-2p)(3p-l) 

(2p-l) P (2-3p) 

l-p 
(l-p)(3p-2) 

2/3-1 



The diffusion coefficient reads 



(1- 


2p)- 


— 2 

P 




(1- 




(2p- 


-1)- 



\<p<\ 
\<P<1 
1<P<1 



0<p<\ 
\<P<\ 

\<p<i 
§<p<i 



(50) 



(51) 



in the symmetric version. 

Generally for the repulsive process with Hamiltonian 



( 48 1 , the current is given by [14] 

p(i-fep) 
i-(fc-i)p 
[(j+i)p-i](i-jp) 

= \ [j-(j+ifp](jP-j+i) 
i-p 

(l-p)(/cp-fc+l) 
(fc-l)p-fe+2 



o< P <l 



k 

< p < i 



(52) 



fc-l 

k 



Here j = 2, 3, . . . , k — 1. The diffusion coefficient charac- 
terizing the models with symmetric hopping is [14] 



D{p) 



[1 - (* - l)p]" 2 
p- 2 

a-p)- 2 

[(fc-l)p-fc + 2]- 2 



0<P<1 
\ <P< 2 

h<P< k T 1 
fc-l 
fc 



(53) 



< p < 1 



In Sees. |III|fV[ we computed the limiting shapes corre- 
sponding to the Ising ferromagnet with NNN interactions 
(k = 2). We now briefly describe limiting shapes arising 
in situations when k > 2. 



2. If a hop would not change the number of NN pairs 
of particles, and it would not increase the number 
of NNN pairs of particles, it is always performed. 

3. Otherwise a hop is never performed. 

This dynamics arises independently of the coupling con- 
stants as long as the inequality Ji > J3 is obeyed. Gen- 
erally the positive coupling constants in ( 48 1 are assumed 
to obey the set of constraints 



Ja > Ja+1 



Jk 



a = 2, 



1 



(49) 



These constraints simplify the analysis since they allow 
us to treat interactions in a lexicographic order. 

The analysis of the repulsive process with Hamiltonian 



(47) is rather cumbersome, yet the results are neat and 
simple. In the totally asymmetric case (only hopping to 



A. Melting of the Quadrant 



We use the same notations as in Sec. IIIII The scaled 
density profile satisfies 



if- 2 /)' + 2(f = 0, 
([l-(fc-l)/]- 2 /')' 



< f < 1 
- Kf = 0, 



< £ < 00 



(54a) 
(54b) 



where I is the scaled position on the one-dimensional lat- 
tice where the density crosses 1/k and the diffusion coef- 
ficient (53) changes its dependence on the density. The 



boundary conditions are 

/(o) = i m 



i /(oo) = (55) 



The same analysis as in Sec. |III| shows that the diagonal 
point on the interface is given by (21 ) with C = — 4/'(0). 



The value of the amplitude depends on k and it can be 
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determined by numerically solving Eqs. (54a|— (54b I sub- 



ject to the boundary conditions (55) 



It is again possible to compute the average area avoid- 
ing a numerical solution of the boundary-value problem, 
Eqs. ( 54a H 54b I and (55), and even without finding £. 



Equations ( 23 ) and ( 25 ) arc still applicable. To compute 
I we multiply Eqs. ( 54a )— ( [54b] ) by £ and integrate the 
former from £ = to £ = £ and the latter from £ = i to 
£ = oo. Using integration by parts we find 



4 / dCC/CC) 

'0 



l/k 



df 



[l-(k-l)f] 



1/2 



l/k 



( 1 
P 



= fc-1 



Thus I = (k — l)/2 and hence the average melted area 
[see Eq. ( 23 )] is given by 



(S t ) = 2(fc - 1)4 



B. Quadrant in a Magnetic Field 



(56) 



Here we consider the quadrant in a magnetic field fa- 
voring the majority phase. Consider first the situation 
when k = 3. Inserting the scaling form ([33]) into Eq. (31 ) 



and using Eq. (50) for the current, we obtain a rather 



simple expression for the scaled density 



1 



R(Z) 



-1/21 



1 [l + (3 + 2Z)" 
\ [l-(3-2Z)- 1 /2] 





Z < -1 
-1 < Z < 
< Z < 1 
Z > 1 



(57) 



For Z < 0, we have p > |(l + 1 



> 



while for 



Z > we have p < |(l - l/VS) < §. Therefore only the 
behaviors of the current in the regions < p < | and 
| < p < 1 matter. 

Generally, the current vanishes at 2k — 1 points, but 
again only the behavior of the current in the regions be- 
tween the first two and last two zeros matter, viz. the 
behavior of J(p) when < p < I and < p < 1. The 
generalization of ( 34 ) and ( 57 ) reads 



R{Z) 



1 

fc-2+[fc+(fc-l)Z]~ 1/2 
fc-1 

l-[fc-(fc-l)Z]~ 1 / 2 
fc-1 





Z < -1 
-1< Z < 
< Z < 1 
Z > 1 



(58) 



Using Eq. ( 35 1 we find that the diagonal point on the 
interface, X^ing 



dZ ■ 



1 



idiag = X k , is given by 
- [k- (k- l)Z]-!/2 



1 



1 



(Vfc + i) 5 



(59) 



Combining Eq. ( 35 1 and ( 58 ) we find the limiting shape 
Y = Y(X) in the range X k < X < 1: 

y = fc-(fc-2)X-2Vfc-l-(fc-2)X 



(fc-2) 2 



In the region < X < X k the limiting shape is found by 
using the mirror symmetry: 



k-(k-2)Y-2y/k-l-(k-2)Y 

x (j^w ( } 

In particular, for k = 3 the limiting shape is given by 

Y = 3-X -2V2-X, (V3 + l)- 2 < X < 1 
X = 3-Y -2V2-Y , 0<X<(V3 + 1)- 2 
This limiting shape is plotted on Fig. |4j 



(62) 



C. Ising Fingers 

In the reference system in which the finger is station- 
ary, the governing equations for the interface are 



v = Vx 




[l-(k- 2)y x ] 

2 



k + 2}- 



0<Vx<k L I 

fc^T < Vx < 1 
1 < y x < k — 1 

k — 1 < y x < oo 



(63) 



Using the same treatment and notation as in Sec. |VJ 
one finds the velocity of the finger 



V 
I 1 



V = 2(k- 1) 



(64) 



The finger is composed of four pieces on the right from 
the tip, < X < e u h < X < £ 2 , h < X < e 3 and 
£3 < X < 1, supplemented by four analogous pieces on 
the left (— 1 < X < 0). The separation points are 



h = 



1 

V' 



1 



ia = 1 



The limiting shape can be expressed through elementary 
functions. The bottom part (0 < \X\ < £{) of the finger 



Y 



\X\ 



k-2 V{k-2f 



\a[l + V{k-2)\X\) (65) 



is close to a parabola, Y ~ ^VX 2 , near the tip. The 
other parts are 

Y = C 1 -V- 1 \n(k~V\X\) (66) 

for l x < \X\ < e 2 , 

Y = C 2 + \X\ + \V{\X\-l 2 ) 2 (67) 

for i 2 < \X\ < £ 3 , and 

Y = C 3 + (k-2)\X\-V- 1 \n(l-\X\) (68) 

for £ 3 < \X\ < 1. The constants Ci, C 2 , C3 are found by 
matching the shapes at common separation points. For 
instance, matching (65) and (66) at \X\ — £1, one finds 



Ci = 



1 



1 



V(k-2) V[k-2) 



2 ln(fc - 1) + - ln(fc - 1) 



1 
V 
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VII. CONCLUDING REMARKS 

We investigated a zero-temperature spin-flip dynam- 
ics of the Ising model on the square lattice with ferro- 
magnetic spin-spin interactions spanning beyond nearest 
neighbors. Specifically, we studied the evolution of a sin- 
gle interface separating two ordered phases. We showed 
that the mapping of an unbounded interface onto a lat- 
tice gas leads to an interesting family of lattice gases 
with repulsive interactions. These repulsion processes are 
tractable [14] , and this allowed us to determine the limit- 
ing shape of the interface in the situation where initially 
the interface is either the boundary of the quadrant, or 
the boundary of a semi-infinite stripe with width greatly 
exceeding the lattice spacing. We found that the limit- 
ing shape of interfaces in the Ising model depend on the 
interaction range. Thus we obtained an infinite family of 
different limiting shapes for the melting quadrant (with 
and without magnetic field) and for the Ising finger. 

An important and feasible extension of our work is to 
consider the shrinking of a droplet of one phase in a sea 
of the opposite phase. This droplet, evolving according 
to a zero-temperature spin-flip dynamics, will eventually 
disappear. In the intermediate regime when the droplet 
is still very large compared to the lattice spacing, but al- 
ready very small compared to the initial size, the droplet 
has an essentially deterministic shape (fluctuations are 
negligible in comparison with the characteristic size of 
the droplet) and the initial condition no longer matters. 
In the case of the classical Ising model, this limiting shape 
has been analytically determined [TT]; see also [3T] and 
the proof of the convergence to the limiting shape [22] . 
It appears possible to generalize these results to include 
long-ranged interactions. 

Another obvious extension is to three dimensions. This 
is a much more challenging domain, e.g., little is known 
about the fate of three-dimensional Ising ferromagnets 
quenched to zero temperature [23] . Limiting shapes 
are also unknown, and even evolution equations govern- 
ing the motion of interfaces haven't been established. 
The mapping of interfaces onto two-dimensional lattice 
gases is possible, but it has not yet led to success. In the 
simplest case of the melting octant in the presence of a 
magnetic field, the governing evolution equation for the 
interface has been proposed in Ref. |24| . It would be in- 
teresting to generalize this conjectural evolution equation 
to the situations with long-ranged interactions. 

We are grateful to S. Redner for discussions and collab- 
oration at an earlier stage of this work. The research of 
J.O. is supported by NSF Grant No. DMR-1205797. 



erning equation is the continuity equation (31). The 



current is given by Eq. (32) in our lattice gas, or by 



J = p(l — p) for the asymmetric exclusion process which 
arises for the standard Ising model with only NN spin 
interactions. The strategy is therefore to solve (3l]) sub- 
ject to the proper initial condition, like Eq. (19) for the 
corner case, and then to use ( 16 1 to extract the limiting 



shape. Here we establish a closed equation describing the 
evolution of the interface. To deduce such an equation 
for the interface y = y(x,t), we re-write (16) as 



y(x,t) = / dzp(z,t) (Al) 

J x—y(x,t) 



Taking the time derivative of (Al), we obtain 



Vt = PVt + j dz p t = py t + J (p) (A2) 
Jx-y 

where on the second step we have integrated pt = —J z - 
Differentiating (Al) with respect to x, we obtain 



Vx = p(y x - 1) 



(A3) 



Hence we can express the lattice gas density through the 
slope of the in terfa ce, p = y x / (y x — 1). Using this result 
we recast Eq. ( A2 ) into a closed evolutionary equation 



Vt = (1 - Vx) J 



Vx 



Vx - 1 



(A4) 



When the current is given by Eq. ( 32 ) , the general equa- 
tion (|A4"1) reduces to Eq. ((38]). 



Without a magnetic field, the governing equation in 
the lattice gas framework is the diffusion equation (JsJ) . 
Taking the time derivative of (Al), we obtain 



Vt = PVt - D(p) 



dp 

dz 



(A5) 



Using dz = dx(l - y x ) and 1 - p = (1- y x ) \ 
Eq. (fA5l) into 



we recast 



Vt = -D(p)p x 



(A6) 



Relation ( A3 1 remains valid and this allows one to reduce 



( A6 ) into a closed evolutionary equation 



Appendix A: Equations Governing the Evolution of 
Interfaces 



Vt = D 



y x -lj (1 - yxf 



(A7) 



Consider first the simpler case of the evolution in a 
magnetic field. In the lattice gas framework, the gov- 



When the diffusion coefficient is given by (12), the gen- 
eral equation AA7h reduces to (|39a|)-(|39b|). 
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